Anisotropic optical conductivities due to spin and orbital orderings 
in LaV0 3 and YV0 3 : First-principles studies 
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The anisotropy of low energy (0~5eV) optical excitations in strongly correlated transition-metal 
oxides is closely related to the spin and orbital orderings. The recent successes of LDA+C/ method 
in describing the magnetic and electronic structures enable us to calculate the optical conductiv- 
ity from first-principles. The LaVC>3 and YVO3, both of which have 3d 2 configuration and have 
various spin and orbital ordered phases at low temperature, show distinct anisotropy in the optical 
spectra. The effects of spin and orbital ordering on the anisotropy are studied in detail based on 
our first-principles calculations. The experimental spectra of both compounds at low temperature 
phases can be qualitatively explained with our calculations, while the studies for the intermediate 
temperature phase of YVO3 suggest the substantial persistence of the low temperature phase at 
elevated temperature. 



> 

O 

O 

-i— > 



T3 

c 

o 
o 



X 



PACS numbers: 75.30.-m, 78.20.Bh, 71.27. +a 



I. INTRODUCTION 



The strong couplings among the spin, lattice and 
charge degrees of freedom in transition-metal oxides 
(TMO) are essentially mediated by the orbital degree 
of freedom (ODF), which plays a crucial role iujiQntrol- 
ling the phases and various physical propertiesoaa. The 
direct observatioax>£orbital structure is difficult, yet sev- 
eral experimentsa'EM have been developed to detect the 
anisotropy induced by spin and orbital orderings. Among 
them, the measurement of anisotropic optical conductiv- 
ityElcl by using polarized light can provide us with useful 
information. The low-energy optical excitations below 
the strong O 2p to transition-metal 3c? charge-transfer 
peak mostly come from the transition-metal 3d-to-3d 
transitions. Those transitions reflect sensitively the spin 
and orbital structures of the system through optical tran- 
sition matrix elements, and show anisotropy when spin 
and orbital are cooperatively orderedQB. Nevertheless, 
the proper analysis of experimental anisotropic optical 
spectra requires the detailed understanding of electronic 
and magnetic structures. 

Both LaVC>3 and YVO3 are insulators with a Pbnm 
orthorhombic unit cell with a fa b ~ c/v2 at room tem- 
perature. There are two 3d electrons per V 3 tsite. Syn- 
chrotron X-ray and neutron diffraction results^ suggested 



that LaV03 undergoes a magnetic phase transition at 
143KEj and a structural phase transition at 140K. The 
low temperature phase has the C-type antiferromagnetic 
(AF) spin configuration (i.e., ferromagnetic (FM) cou- 
pling along the c-axis and AF coupling in the afo-plane). 
Due to the Jahn- Teller distortion, one of the V-0 bond 
is longer than other two in the VOg octahedron. The 
longer bond lies in the a6-plane and its orientation al- 
ternates for the neighboring V sites in the plane. The 
stacking of the longer bonds also alternates along the 
axis. Hereafter we call this kind of distortion the G-type 
Jahn- Teller distortion. This low temperature phase has 
P2i/a crystal symmetry. 

The structural . and., magnetic phases for YVO3 are 
quite complicatecOEJiiil With lowering the tempera- 
ture, this compound first undergoes a structural phase 
transition at 200K from a disordered phase to the G-type 
Jahn- Teller distorted structure, the same structure as the 
low temperature phase of LaV03. Then the C-type AF 
ordering develops at 116 K (Tjvi). With further lowering 
the temperature to 77K (T/V2), another structural and 
magnetic phase appears. The low temperature (<77K) 
phase recovers the Pbnm crystal symmetry with the C- 
type Jahn- Teller distortion (i.e., the longer bonds stack 
along the c-axis in the same orientation rather than alter- 
nately) . The magnetic structure of this low temperature 
phase is G-type AF structure (i.e., AF coupling both in 
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the a&-plane and along the c-axis). Strong temperature- 
induced ^Magnetization reversal can be observed at T/vi 
and T N2 £-3. 

Sawada and Terakura studied the electronic and mag- 
netic structures of LaVC>3 and YVO3 by using the full- 
potential linearized augmented-plane-wave (FLAPW) 
method. They first compared the resultsEJ between LDA 
(local density approximation) and GGA (generalized gra- 
dient approximationiy, and then applied the LDA+J7 
methodEJ to LaV0 3 B. They found that the LDA and 
GGA are not sufficient, and LDA+C7 is necessary in or- 
der to predict the low temperature phase correctly. Re- 
cently, the polarized optical conductivities otLaVOs and 
YVO3 were measured by Miyasaka et al.lla and very 
clear anisotropy was observed. To understand the im- 
plications of those spectra, we developed the plane-wave 
pseudopotential method based on LDA+C7 to calculate 
the optical conductivity and studied the effects of spin 
and orbital orderings on the anisotropy. The experimen- 
tal spectra, of both LaVO;; and YVO3 at low tempera- 
ture can be qualitatively explained by our calculations, 
while the studies for the intermediate temperature phase 
(77K < T < 116K) for YV0 3 suggest the complication 
of magnetic and orbital structures. In Section II of this 
paper, we will describe our method, and the results are 
discussed in Section III. 



II. METHOD 

The Vanderbilt type ultra-soft pseudopotentialEl is 
useful not only for efficient calculations for transition- 
metal oxides but also for implementing the LDA+Z7 
method to treat effects of strong correlation. In the 
LDA+J7 method, the strong Coulomb interaction is ex- 
plicitly taken into account in the subspace of localized 
orbitals through a Hartcree-like scheme. The detailed 
description of our LDA+C/ scheme was given in Ref. 0. 

The inter-band optical conductivity is calculated from 
the converged Kohn-Sham wave functions |V>nk) aad. 
eigen values E n (k) by using the following Kubo formulae^ 
(in Ry units): 
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where a and [3 {=x, y, z) are indices for directions, w 
is the excitation energy, V is the volume of the unit 
cell, n and m are band indices, f n k is the Fermi dis- 
tribution function, u! mn — E m (k) — E n {k) and 8 is the 
lifetime broadening (<5=0.01Ry in this work), 7r" m = 
(ipnk\(— iV Q )|?/'mk) are the matrix elements of the mo- 
mentum operator. The calculations for the matrix ele- 
ments 7r" m require the all electron wave-functions ip n \i, 



which can ha-obtained from the following core compen- 
sation formEilE3: 



|^nk) = |0„k) + ~ l&» W«0 (2) 



where 0„k are the pseudo-wave-functions obtained from 
the self-consistent pseudopotential calculations, i is the 
index for atomic orbitals, ipi and <f>i are atomic all- 
electron and pseudo- wave functions respectively, (3i are 
the localized functions as defined in Ref.|ljJ. Therefore, 
the matrpi .elements of momentum operator can be ob- 
tained aaSEj: 
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Practically, the core contribution (second term of 
eq.(3)) can be calculated in the pseudopotential gener- 
ation process and stored as input of self-consistent cal- 
culations. The present processJias been well checked by 
using LaMn03 as an exampleLM Throughout the cal- 
culations, 30Ry has been used for the cutoff energy of 
plane-wave expansion, and we use (6x6x4) mesh for 
the k-points in the linear tetrahedron method with the 
curvature correction. The parameter U e ff in LDA+C7 
scheme is chosen to be about SiJ^-ip. order to repro- 
duce the experimental band gaps 



III. RESULTS AND DISCUSSIONS 
A. Electronic structures 

The unit cell and coordinates are defined in Fig. 1 . Four 
V atoms, say V-l, V-2, V-3 and V-4, are included in the 
unit cell. Following the common convention, we defined 
the x, y and z directions as the [110], [110] and [001] 
directions of the unit cell respectively. Three structures, 
LaV0 3 at 10KB, YVO3 at 65K1I2, and YV0 3 at 100KEI 
are treated here, which have the G-type, C-type and G- 
type Jahn- Teller distortions, respectively. For each fixed 
structure, relative stability is studied among four kinds 
of collinear magnetic structures: C-type AF, G-type AF, 
A-type AF(i.e., FM layers coupled antiferromagneticaly 
along the c-axis), and FM spin ordering (SO) states. 

Each of the above Jahn- Teller distortions is accom- 
panied by the same type of orbital ordering (OO) and 
stabilizes a specific related magnetic structure. In the 
present case, our calculations suggest that the G-type 
and C-type OO favor the C-type and G-type AF SO, 
respectively. This result is consistent with experim ental 
observations and also with the previous-calculationsLJtll 
and unrestricted Hartree-Fock studiesc.1 Now we con- 
centrate on the discussion for the electronic structures of 
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FIG. 1: The unit cell and coordinates used in the calculations. 
Four V sites (black spheres), labeled as V-l, V-2, V-3 and V-4, 
are included in the unit cell. 



4j,4*; v " 3: 4 y 4 z ; V-4: d\ y d\ x . As for the 00 and SO 
in the a6-plane, the present situation is the same as those 
of A) and C). However, the 00 and SO along the c-axis 
are different. In the present case, all the V atoms along 
a given c-axis have the same orbital occupation, while 
those belonging to the nearest neighbor c-axes have dif- 
ferent orbital occupations. For instance, the d yz orbital is 
occupied for both V-l and V-3 sites. This C- 00 leads to 
stability of the G-SO through AF superexchange. Com- 
paring LaV03 with YVO3, we can find strong mixture of 
orbitals in YVO3. This is due to the enhanced structural 
distortion in YVO3 compared with LaV03. The V-O-V 
angles are about 157 degrees in LaV03, while they are 
reduced to about 144 degrees in YVO3 . 

Similar analysis was made for all other magnetic phases 
of these three structures. Among several results ob- 
tained, we noticed the following interesting and possibly 
important feature: for a fixed crystal structure, a change 
in the magnetic structure will not cause any significant 
change in the orbital ordering. It is likely that the orbital 
structure is mostly determined by the specific structure 
distortion in the present case. 



B. Optical conductivities 



the lowest energy magnetic state of each structure. They 
are A) LaV0 3 10K with G-00 and C-SO; B) YV0 3 65K 
with C-00 and G-SO; C) YV0 3 100K with G-00 and 
C-SO. The obtained magnetic moments for these three 
states are all about 1.7/is/V. The calculated electronic 
projected densities of states (PDOS) are summarized in 
Fig. 2, where different orbitals are shown by different col- 
ors. The 0-2p and V-e 9 states are not plotted. 

For A), the occupied electronic configuration can be 
nominally expressed as V-l: 4y4z' V-2: d xy dl x ; V-3: 

dl y dl x ; V-4: d x d^ z . ^ )ne °^ tw0 e l ec t rons on eacn 
V site will occupy the d xy orbital, and another electron 
will occupy the d yz or d zx orbital alternately from one V 
atom to its neighboring V atoms, resulting in G-00. The 
C-SO is manifested as the up spin state for V-l and V-3, 
and the down spin state for V-2 and V-4. Due to devia- 
tion from the cubic structure, certain degree of mixture 
of orbitals can be observed: for example, the small oc- 
cupation of d\ x orbital (green line) for V-l site. Note 
here that, in our calculation, for each fixed structure 
and spin configuration the orbital polarization is fully 
relaxed self-consistently. The calculated electronic and 
magnetic structures of YVO3 in the intermediate tem- 
perature (100K) phase (panel C in Fig. 2) can be under- 
stood in the same way as the low temperature phase of 
LaV03 in terms of the same spin and orbital orderings. 
This can be seen from the similarity of the panel A) to 
the panel C) in Fig. 2. 

For the low temperature phase of YVO3 (panel B in 
Fig. 2), the situation is different. The nominal electronic 
occupations can be understood as V-l: d x d yz ; V-2: 



The results for optical conductivity are summarized 
in Fig. 3. Experimentally, two peak structures (a and 
(3) are observed for E//c in both compounds, with peak 
(3 located at higher energy than peak a. For LaV03, 
peak a has a very sharp structure at low temperature 
and is significantly suppressed above the transition tem- 
perature. Peak (3 is quite weak and almost temperature 
independent. In YVO3, on the other hand, peak (3 has 
a larger weight with strong temperature dependence and 
peak a is weak with only little temperature dependence. 
For both LaV03 and YVO3, the E//a spectra are quite 
broad and almost temperature insensitive. 



1. LaVOs low temperature phase: 

For LaV0 3 at 10K, the G-00 with C-SO state is the 
most stable state. We first analyze the case of E//c, 
where we consider the inter-site transition between V-l 
and V-3 as the transition dipole is along the c-axis. The 
calculated optical conductivity for C-AF in Fig. 3 shows 
clearly a sharp peak corresponding to the experimental 
peak a. Looking back to the calculated PDOS in Fig. 2, 
and taking the transitions from V-l site to other sites as 
examples, we can assign the sharp peak a for E//c to the 
transition from V-l-c£L to V-3-4, as denoted by the solid 
red arrow in Fig. 2. This transition is not allowed if the 
magnetic structure along c-axis is AF, like in the G-AF 
and A-AF states, where no peak is obtained around the 
same position of peak a. Therefore, the peak a should be 
suppressed above the magnetic transition temperature as 
observed in experiment. On the other hand, peak (3 seen 
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A) LaVO, 10K B) YVCX 65K C) YVO, 100K 
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FIG. 2: The calculated projected densities of states (PDOS) for the lowest energy states of three systems. They are A) 
LaVOs 10K structure with G-OO and C-SO; B) YVO3 65K structure with C-OO and G-SO; C) YV0 3 100K structure with 
G-00 and C-SO. The lines with arrow indicate the optical transition paths from V-l site to other sites. All the PDOS and 
paths for different orbitals are indicated by different colors. The solid arrows show the real transition in this spin and orbital 
configuration, while the dashed arrows show the possible transition in different magnetic structures. The positive and negative 
PDOS values mean t-spin and J.-spin respectively. See the text for detailed explanations. 



in experiment does not exist in the calculated spectra for 
C-AF. Since the spectral weight of peak f3 is very small 
and temperature insensitive in experiment, we speculate 
that this peak may come from certain component of im- 
perfect spin ordering along the c-axis. For example, let 
us reverse the spin of V-3 site in Fig.2A. The transition 
from V-l-d) yz to V-3-(iJ z located at higher energy (de- 
noted by the red dashed arrows in Fig. 2) will be allowed 
and contribute to a sharp peak near the position of peak 
(3, as shown in the spectra for the G-AF and A-AF state 
for E//c. Before moving to YVO3, we briefly discuss the 
case of E//a, where the intersite electronic transitions 
between V-l and V-2 are important. The two transitions 
from V-l-dl to Y-2-d\ (denoted as solid blue and 
red arrows in Fig. 2) are allowed for the a6-plane spec- 
tra, which show a broad structure due to the broad band 
width of unoccupied minority spin states. These a6-plane 
transitions produce a structure in E//a curve at the po- 
sition of peak (3 of E//c curve. Due to the doubling of 
transition paths, one may expect a strong a6-plane spec- 
tral weight. However, by calculating the matrix elements, 
the actual a6-plane spectrum is not so much weighted due 



to the longer bonds in the a6-plane. The shoulder in the 
experimental E//a curve at about 1.5 eV, i.e., the en- 
ergy of peak a in the E//c curve may originate from the 
spin disordering in the a6-plane. First we note that a 
very tiny peak is seen at this energy in the C-AF ground 
state configuration. This structure mainly comes from 
the transition from the occupied V-l-d]^ component to 
a small N-2-d\ x component just above the Fermi level. 
However, this structure is too weak to explain the shoul- 
der at 1.5 eV. If the spin state between V-l and V-2 may 
become parallel like in A-AF and FM, the transition cor- 
responding to the red dashed arrow may become allowed 
and produces a peak as shown in Fig. 3. 

From the above assignment, some important parame- 
ters can be estimated by using the Hartree-Fock model. 
The peak position of a should corresponds to energy 
U — J, and peak (3 corresponds to U + J, where U and J 
are on site Coulomb and exchange parameters. The esti- 
mated values for U and J are 2.2eV and 0.6eV from our 
calculations, and about 2.1eV and 0.32eV from experi- 
ments for LaV03. The smaller J parameter estimated 
from experiments may suggests the spin canting or flue- 
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A) LaV0 3 10K B) YV0 3 65K q YV0 3 100K 
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FIG. 3: The calculated optical conductivities for various phases. Three structures, A) LaV03 at 10K, B) YVO3 at 65K, and C) 
YVO3 at 100K, are calculated. For each fixed structure, results for four different magnetic configurations are plotted together 
with experimental results. The figure enclosed by the red border line corresponds tei the most stable spin structure. The 
experimental curves in A), B) and C) are taken at 10K, 67K, and 110K respectively^. The solid lines are results for E//c, 
while dashed lines for E/ /a. See the text for the detailed explanations. 



tuation at elevated temperature, especially for YVO3 as 
discussed below. 



2. YVO3 at low temperature (T<77K): 

Now let us discuss the low temperature phase of YVO3. 
In this case, the C-00 with G-SO state is the most sta- 
ble state. As the results, the strong peak (3 for E//c 
can be explained as the transition from V-l-d yz to V-3- 
d yz as noted by the solid red arrow in Fig.2B with energy 
U + J. This transition is not allowed for the FM coupling 
along the c-axis as in the case of C-AF or FM state, and 
should show strong temperature dependence as observed 
in experiment. As for the small peak a, the situation for 
C-00 is different from the case of LaV03 with G-OO. 
In the present case of C-OO, the significant mixture of 
different orbitals under strongly distorted V06 environ- 
ment is crucially important to produce this small peak a 
located at energy U — J (corresponding to the transition 
as noted by the green arrow in Fig.2B). Let us explain 



the situation in more details. If there is no orbital mix- 
ing (the inter-site hybridization for the same orbitals still 
exists), d yz is fully occupied and d zx is empty in the ma- 
jority spin state both at V-l and V-3. The unoccupied 
d\ x at V-l hybridizes with unoccupied d\ x (not d yz ) at 
V-3 to produce a small d\ x component at V-3 site at 
the energy of majority spin empty state. In this case, 
the weight of the transition from Y-\-d\ x to V-3-d| x is 
zero because of the zero occupation of \-\-d\ x state, and 
the small peak a should not exist. However, as already 
mentioned the orbital mixing between d yz and d zx is sig- 
nificant due to the lattice distortion and produces the 
partial occupation of Y-l-d\ x state (green line in Fig.2B). 
This will contribute to the small peak a. We should em- 
phasize here that by keeping the C-00 any magnetic 
coupling along the c-axis will produce a peak around a 
only through the orbital mixing and therefore the peak 
a will never be strong, in contrast to the case of LaV03. 
Nevertheless, those states with different orbital ordering, 
like the intermediate temperature phase with G-00 and 
C-SO, can contribute to the peak structure around a as 
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shown in Fig.2C. 

3. YV0 3 at intermediate temperature (77K<T<116K): 

The observed spectra for the intermediate phase of 
YVO3, which has G-00 and C-SO, is hard to be ex- 
plained by our calculations. Our calculations suggest a 
picture similar to that of the low temperature phase of 
LaVC>3, because the two systems have the same spin and 
orbital orderings. However, the observed spectra is quite 
similar to the low temperature phase of YVO3, which has 
different spin and orbital orderings. In particular, the ob- 
served peak (3 is stronger than peak a and temperature 
dependent, while the calculations suggest the opposite. 
One possible scenario to explain the observed picture is to 
suggest the persistence of low temperature phase (C-00 
and G-SO) at elevated temperature, or at least certain 
mixture of other magnetic states, like G-00 with G-SO 
or G-00 with A-SO states. This argument can be sup- 
ported by the following facts. First, the small J param- 
eter (about 0.25eV) estimated from experiment suggests 
the possibility of strong canting or magnetic fluctuation 
at elevated temperature (100K). Second, the recent neu- 
tron experiments on this compound actually suggests 
the existence of G-type spin canting and significant re- 



duction of magnetic moment at 100K. Nevertheless, ac- 
cording to our calculations, to be able to reproduce the 
observed shape, substantial (more than 50 percent) mix- 
ture of low temperature phase is required. 

In conclusion, we developed the first-principles plane- 
wave-pseudopotential method based on LDA+J7 to cal- 
culate the inter-band optical conductivity. We calcu- 
lated the anisotropic optical conductivities of LaV03 
and YVO3 for different orbital and spin ordered phases, 
and studied the effects of spin and orbital ordering on 
the anisotropy. The experimentally observed spectra for 
both LaV03 and YVO3 at low temperature can be qual- 
itatively explained by our calculations. On the other 
hand, our calculation for the intermediate temperature 
phase for YVO3 suggests that substantial mixture of low 
temperature phase may persist at elevated temperature. 
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